Color superconductivity vs. pseudoscalar condensation in a three-flavor NJL model 



Harmen J. Warringa and Daniel Boer 
Department of Physics and Astronomy, Vrije Universiteit, 
De Boelelaan 1081, 1081 HV Amsterdam, The Netherlands 



Nordita, 



Jens O. Andersen 
amsvej 17, DK-2100 Copenhagen 0, Denmark 
(Dated: February 2, 2008) 



in 
o 
o 

(N 



ON 



> 

i> ■ 

^t- : 
o . 
in ■ 

o ; 

^ : 

a: 
i 

<D . 
43 ■ 



X 



We calculate numerically the phase diagram of the three-flavor Nambu-Jona-Lasinio model at 
zero and finite temperature as a function of the up, down, and strange quark chemical potentials. 
We focus on the competition between pseudoscalar condensation and color superconductivity. We 
find that the two types of phases are separated by first-order transitions. 



I. INTRODUCTION 

In recent years a huge effort has been made in order 
to understand the phase diagram of strongly interacting 
matter, both experimentally as well as theoretically. Col- 
lisions at the Relativistic Heavy-Ion Collider (RHIC) at 
BNL and the Large Hadron Collider (LHC) at CERN al- 
low the experimental study of hadronic matter at energy 
densities exceeding that required to create a quark-gluon 
plasma. The energy and baryon densities of these exper- 
iments correspond to temperatures up to 200 MeV and 
a baryon chemical potential in the range 0-600 MeV. 

At zero baryon chemical potential, lattice calcula- 
tions suggest the existence of a transition from ordinary 
nuclear matter to a quark-gluon plasma at temperatures 
of around 160 MeV (see for example Ref. pj). While 
there has been progress in generalizing lattice calcula- 
tions to finite baryon chemical potential |2J,y,|^,|5|], highly 
nontrivial problems remain to be solved, and in practice 
one is restricted to small baryon chemical potentials. 

Since one cannot apply lattice gauge theories to the 
region of the phase diagram where the baryon chem- 
ical potential is large, one must employ other meth- 
ods. One possibility is perturbative QCD. Asymptotic 
freedom guarantees that this method can be applied at 
asymptotically large chemical potentials. If the densi- 
ties are too small to use perturbative QCD, one may 
apply effective theories. Models such as the instanton 
liquid model, random matrix models, and Nambu-Jona- 
Lasinio (NJL) models have been applied to the study of 
the QCD phase diagram at finite temperature and densi- 
ties. Despite their shortcomings, it is expected that such 
models do describe the qualitative features of the QCD 
phase diagram, in regions not accessible to perturbative 
or lattice QCD. For example, the existence of a criti- 
cal point at finite baryon chemical potential has been 
predicted using effective models [HQ- Furthermore, it 
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is expected that quark matter at high baryon chemical 
potentials and low temperatures can be in a two-flavor 
color-superconducting pha se ( 2SC) d, |j| ^| or a color- 
flavor locked (CFL) phase [TlLIT^. In these phases gaps 
on the order of 100 MeV arise. 

Many results have been obtained for equal up, 
down and strange quark chemical potentials. However, 
this may not be directly relevant for heavy-ion colli- 
sions or compact stars. For example, to enforce elec- 
tric and color neutrality and weak equilibrium in com- 
pact stars different flavor and color chemical potentials 
have been introduced 0, This gives rise to a 

more complicated phase diagram in which one also finds 
new 2SC-lik e Il3 and s apless superconducting phases 
[H [H H ElM ETEi E|. Similarly, in heavy- 
ion collisions, a difference between the quark chemical 
potentials arises if the number densities of the different 
quark flavors are not the same. This difference can cause 
interes ting observable effects such as two critical end- 
points |24l l25j|. However, instanton induced interactions 
tend to suppress this effect [2f|. 

In addition to a more complicated structure of the 
superconducting phases, different chemical pot entials can 
also trigger pseudoscalar condensation [23, IH, |2!j • This 
has been confirmed on the lattice at zero baryon but finite 
isospin chemical potential [3(J • Depending on the flavors 
involved, the charged pion, neutral or charged kaon fields 
may acquire a vacuum expectation value. Pseudoscalar 
condensation in the two-flavor NJL model has been stud- 
ied in Refs. [25ll3l| . as a function of the different chemical 
potentials at zero and finite temperature. An extension 
to three flavors was carried out in Ref. [13 • 

In this paper, we discuss the phase diagram of the 
three-flavor NJL model as a function of the different 
chemical potentials including both pseudoscalar conden- 
sation and color superconductivity. We did not apply 
electric or color neutrality conditions. At zero tempera- 
ture, pseudoscalar condensation is p ossible if \fi u — > 
m w [23 or if \fi Ul d — Ms I > m K \2M ■ On the other hand, 
color superconducting phases occur if the chemical poten- 
tials are large and approximately equal. Therefore, one 
can imagine scenarios where for example \i u ~ [id (the 



Fermi surfaces of the u and d quark should be sufficiently 
close for Cooper pairing to occur) and /i u w — /i s (the 
Fermi surfaces of the u and s should be sufficiently close 
for kaon condensation to occur), with \/j, u ,d ~ Ms I > m-K- 
In such a case a 2SC phase is competing against a phase 
in which kaons condense. From our calculations it fol- 
lows that a coexistence phase of pseudoscalar condensa- 
tion and color superconductivity does not occur for the 
parameters chosen and that these phases are separated 
by a first-order transition. However, we do not exclude 
that other choices of parameters may lead to such a coex- 
istence phase, just as coexistence of color superconduc- 
tivity and chiral symmetry breaking may occur in the 
NJL model for specific ranges of parameters Here 
we refer to a coexistence phase to mean a phase in which 
two condensates are nonzero simultaneously. 

In this work we study pseudoscalar condensation in 
the quark anti-quark channel. We have not taken into 
account the pseudoscalar diquark interaction. This in- 
teraction is suppressed relative to the scalar diquark in- 
teraction due to instantons 0. However in absence of 
instanton interactions, if one neutralizes the bulk matter 
with respect to color and electric charges it is possible to 
have pseudoscalar diquark condensation with rather large 
gaps j2J] ■ Pseudoscalar diquark condensation in the NJL 
model is similar j^H to pseudoscalar condensation in the 
CFL phase studied with effective chiral models in Refs. 
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We emphasize that in this paper we do not impose 
electric or color neutrality conditions, but that this would 
qualitatively affect the phase structure, leading for ex- 
ample to the observation that in a macroscopic volume 
of quark matter the 2SC phase is energetically disfa- 
vored 0, 0|. Therefore, our results do not address 
the qualitative features of the QCD phase structure with 
neutrality constraints imposed. 

The paper is organized as follows. In Sec. II, we 
briefly describe the NJL model and the choice of param- 
eters. In Sec. Ill, we discuss the calculations and some of 
its technical aspects. In Sec. IV, we present our results, 
and in Sec. V we conclude. 



pressed the color, flavor, and Dirac indices of the fermion 
fields ip for notational simplicity. The mass matrix Mq is 
diagonal and contains the bare quark masses tjiq u , mod 
and ihqs- The matrix \x is also diagonal and contains the 
quark chemical potentials fi u , fid and fi s . We use the met- 
ric = diag(H ) and the standard representation 

for the 7-matrices. The quark-antiquark interaction part 
of the Lagrangian density is 



Cr, 



G 



(ipX a ip) 2 + (ipXai'ysipy 



(2) 



The matrices A a are the 9 generators of U(3) and act 
in flavor space. They are normalized as TrA Q A;, = 2S a b- 
The diquark interaction term of the Lagrangian density 
is given by 



(3) 



where A, B £ {2,5,7} since only the interaction in the 
color and flavor antisymmetric triplet channel is attrac- 
tive. The matrices t a are the generators of U(3) and act 
in color space. Their normalization is Tr t a tb = 28 a b- To 
remind the reader, the antisymmetric flavor matrices A2, 
A5 and A7 couple up to down, up to strange and down to 
strange quarks, respectively. The charge conjugate of a 
field ip is denoted by tp c — C^i T where C — 17072- The 
coupling strength 3G/4 of the diquark interaction is fixed 
by the Fierz transformation. However, some authors dis- 
cuss the NJL model with a different diquark coupling 
constant (see for example Ref. for a comparison). 

The results that will be presented below are obtained 
with the following choice of parameters 

mou = mod = 5.5 MeV , mo s = 112 MeV 

G = 2.319/A 2 , A = 602.3 MeV. (4) 

This choice of parameters gives rise to constituent quark 
masses M u = M d = 368 MeV and M s = 550 MeV 



III. THE EFFECTIVE POTENTIAL 



II. THE NJL MODEL 

In the NJL model [4^, one treats the interaction 
between the quarks as a point-like quark color current- 
current interaction. This naive approximation to QCD 
works very well in explaining various low-energy observ- 
ables such as hadron masses ^3|. By applying several 
Fierz transformations to the current-current interaction 
(see for example Rcf. 44] ) and including only terms 
which give rise to attractive qq and qq channels, one ob- 
tains the following Lagrangian density 

C = $ {vfdp - M + a*7o) 1> + Cqq + C qq , (1) 

where the quark-antiquark term, Cq q , and the diquark 
interaction term, C qq , are defined below. We have sup- 



To obtain the phase diagram of the NJL model, we 
first introduce 18 real condensates a a and f3 a , and 9 com- 
plex condensates Aab as follows 

a a = -2G($\ a iP) , (5) 
/3 a = -2G <^A a z 75 ^) , (6) 

A AB = ^G(tp T t A \ B C l5 ^) . (7) 

We will assume that all condensates are space-time in- 
dependent. The crystalline Larkin-Ovchinnikov-Fulde- 
Ferrell (LOFF) phase 0] will not be considered here. 
The next step is to apply a Hubbard-Stratonovich trans- 
formation to eliminate the four-point quark interactions 
to make the Lagrangian quadratic in the quark fields. 
Introducing a two-component Nambu-Gorkov field ^ T = 



2 



(ip T , ip^) I V% allows for straightforward integration over 
the quark fields. After going to imaginary time, the ther- 
mal effective potential V in the mean-field approximation 
reads 



V 



4G 



3G 



p =(2n+l)7rT 

where K is a 72 x 72 matrix 



- V / ^logdetX, (8) 



if 



and 



-A* AB t A O As ® 75 



(9) 



Z?i = 1/(8 (i7op + 7ipj) - (g) 70 

-(Mo + aaAa)<g>ld-/?aA Q ®i75 , (10) 

2? 2 = 1/ ® (*7oPo + 7»Pi) + A* ® 7o 

-(Af + a a A^) ® l d - /3 a Aj <8 «75 ■ (H) 

The matrix 1 is the identity matrix in color (c), flavor 
(/), or Dirac (d) space. 

The values of the condensates and the phase diagram 
are determined by minimizing the effective potential V 
with respect to the condensates. To make the minimiza- 
tion procedure easier, one can take advantage of the fact 
that certain condensates must vanish. Firstly, applica- 
tion of QCD inequalities [2?], shows that in QCD at 
zero chemical potential chiral symmetry breaking can- 
not be driven by parity-violating condensates of the type 
("0175^). Outside the phase in which diquarks condense, 
we found numerically that this is also correct in the NJL 
model at finite chemical potentials. Therefore, 0o, 03 
and 0s are put to zero. Secondly, although perturbative 
one-gluon exchange cannot distinguish between 0k and 
at condensation with k S {1,2,4,5,6,7}, pseudoscalar 
condensation is favored due to the instanton interaction 
|27j . We therefore set all a fc 's with k <E {1,2,4,5,6,7} 
to zero. We found numerically that this is correct, de- 
spite the fact that the model we consider does not include 
instanton interactions. 

One can further simplify the minimization procedure 
by using the symmetries of the NJL model. In absence 
of quark masses and chemical potentials, the Lagrangian 
density has a global SU(3) C x U(3) y x U(3) A symmetry. 
Due to the non-vanishing quark masses, the symmetry 
is broken down to SU(3) C x U(3) y . Since we consider 
different quark masses and finite chemical potentials, the 
symmetry of the Lagrangian density is further reduced 
to SU(3) C x U(l) u x U(l) d x U(l) s . The vacuum man- 
ifold is invariant under the same transformations as the 
Lagrangian density, so applying a U(l)-flavor transforma- 
tion to all condensates leaves the free energy invariant. 
Therefore, using the U(l)-flavor transformations one can 
choose the pseudoscalars to condense in the 02, 05, and 
01 channels, and set 0\, 04 and 0e to zero. The phase in 



which 02, 05 or 0-j is non- vanishing is called the ir + /tt~ , 
K°/K° or K + /K~ condensed phase, respectively. 

Because of the global SU(3) C symmetry, one can also 
rotate away several diquark condensates. Without loss 
of generality, we can minimize with respect to A22, A25, 
A55, A27, A57 and A77. In principle, all six diquark con- 
densates can have a phase. It is always possible to remove 
two of them by using the two diagonal SU(3) C transfor- 
mations. As long as there is no pseudoscalar conden- 
sation, one can use the U(l)-flavor symmetries to rotate 
away three other phases. As a result either A25, A 55 , A27, 
or A57 has a phase 0|. However, this reduction is not 
completely possible if pseudoscalar condensation occurs. 
By choosing the pseudoscalars to condense in the 02 , 05 
and 07 channels, one breaks the U(l)-flavor symmetry. 
Hence if pseudoscalar condensation arises in one chan- 
nel, one can in general rotate away one phase less in 
the diquark sector. If it occurs in more channels, two 
phases less can be rotated away. However, numerically 
we find that allowing for a complex phase leads to di- 
quark condensation only in the A22, A55 and the A77 
channels. The A25, A27, and A57 diquark condensates 
do not arise or can be rotated away. Moreover, we find 
that pseudoscalar condensation in the quark-antiquark 
channel does not coexist with color superconductivity, 
such that one can always take the diquark condensates 
to be real. 

The different possible color-superconducting phases 
are named as follows |23| 



A 22 


^0, 


A 55 


^ 0, A77 


^0 


CFL , 


A77 


= 0, 


A 22 


+ 0, A 55 


^0 


uSC , 


A 55 


= 0, 


A 22 


£ 0, A77 


^0 


dSC , 


A 22 


= 0, 


A 55 


+ 0, A77 


^0 


sSC , 


A 22 




A 55 


= 0, A77 
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2SC , 


A 55 




A 22 


= 0, A77 


= 


2SCus 


A77 


y^0, 


A 22 


= 0, A 55 


= 


2SCds 



(12) 



To calculate the effective potential in an efficient 
way, one can multiply the matrix K with diag(l c ® 1/ £g> 
70, l c ® 1/ ® 7o) which leaves the determinant invariant. 
In this way, one obtains a new matrix K' with ip$s on 
the diagonal. By determining the eigenvalues of the ma- 
trix K' with pq — 0, one can reconstruct the determinant 
of K for all values of po which is namely rii=i (Ai ~ Wo)- 
After summing over Matsubara frequencies, one finds 



T logdet^ = 

p = (2n+l)irT 



72 

E 

i=l 



y+Tlog 



-Ai /T 



(13) 

All that remains in order to determine the effective po- 
tential is to integrate over three-momentum p up to an 
ultraviolet cutoff A. 

The speed of the calculation of the effective poten- 
tial depends heavily on how fast one can compute the 
eigenvalues. There are several ways to speed up the cal- 
culation. Firstly, the determinant of K does not depend 



3 



on the direction of p. Therefore, one can choose p to 
lie in the z-direction. Together with the choice of the 
non- vanishing condensates mentioned above, this implies 
that K'(po — 0) becomes a real symmetric matrix, which 
simplifies the calculation of the eigenvalues. Secondly, 
one can interchange rows and columns of K' without 
changing its determinant. By doing so, one can bring 
K' in a block-diagonal form. One can then determine 
the eigenvalues of the blocks separately which is signifi- 
cantly faster since the time needed to compute eigenval- 
ues numerically scales cubically with the dimension of the 
matrix. In the most general case with diquark condensa- 
tion, one can always reduce the problem to two 36 x 36 
matrices. Moreover, if there is no diquark condensation, 
but only pseudoscalar condensation, the problem can be 
further reduced to computing the eigenvalues of two 6x6 
matrices. 

We determined the eigenvalues using LAPACK rou- 
tines |47| . After numerical integration over three- 
momentum p up to the cutoff, the condensates were de- 
termined by minimizing the effective potential using MI- 
NUIT |48j . To be certain that the minimization pro- 
cedure did not end up in a local minimum, we always 
checked the continuity of the minimized effective poten- 
tial as a function of chemical potentials and/or temper- 
ature. 



IV. PHASE DIAGRAMS 

In this section, we present our results for the phase 
diagrams of the NJL model with u, d, and s quarks. We 
plot the phase diagrams as a function of the chemical 
potentials and temperature. To determine the locations 
of the phase boundaries, we examine the behavior of the 
condensates. If a condensate jumps discontinuously the 
transition is first-order, and this is indicated by a solid 
line. If its derivative has a discontinuity, the transition 
is second order, and this is indicated by a dotted line. 
If a condensate changes rapidly in a narrow range with- 
out vanishing there is a smooth cross-over, and this is 
indicated by a dashed-dotted line at the point were the 
condensate varies maximally. 



A. fi 3 = and T = 

In Fig. n we display the phase diagram of the NJL 
model for fi s = and T = as a function of fi u and fid- 
Outside the 2SC phase (q), our results agree qualitatively 
with the two-flavor calculations of Ref. [3l| (see their Fig. 
1), where color-superconducting phases were not taken 
into account. Moreover, in Ref. |3jJ different parameters 
were used, in addition to a form factor which mimics 
asymptotic freedom. 

One can clearly see that the phase diagram is sym- 
metric under reflection in the origin. This is because 
the free energy is invariant under the transformation 




liu (MeV) 

FIG. 1: Phase diagram for fi a — and T = as a function of 
fi u and fid- First and second-order transitions are indicated 
by solid and dotted lines, respectively. The letters denote the 
different phases, where a: uu + dd + ss, c: uu + ss, d: dd + 
ss, g: ss, i: 7v + /tt~ + ss and q: 2SC + ss. 



(Mui A*d, /■*«) -> {-Hu,-Hd,-Hs) horn the symmetry be- 
tween particles and anti-particles. Fig.^is also symmet- 
ric under interchange of u and d, because of the choice of 
equal up and down quark masses. This gives rise to the 
symmetry of the phase boundaries with respect to the 
diagonals. 

In general, horizontal and vertical lines in the phase 
diagrams arise if the pairing of one type of quark is not 
changed after a transition. In this case, the location of 
the phase boundary is determined by the properties of 
other quarks. Therefore, changing the chemical poten- 
tial of the unchanged quark species cannot have a big 
influence on the location of the phase boundary. This 
results in the horizontal and vertical lines. For T = 0, 
one always finds these lines near the values of the con- 
stituent quark masses, i.e. fi u « M u , fid « Md and 
pb s w M s (see for example Ref. 0]). The diagonal lines 
arise because at T — pion condensation can occur if 
\fjL u - fi d \ > TO* = 138 MeV H3. 

The diagram shows that if the chemical potentials 
are different, the transition to the color-superconducting 
phase (q) remains first order as was concluded in Ref. 
|49j . Moreover, one can see from Fig. that if fi u ^ fid 
it is possible to go through two first-order transitions be- 
fore entering the 2SC phase (q) (similar to the situation 
discussed in Ref. |25| without color superconductivity). 
We observe that to have such a scenario at zero tem- 
perature, a minimum difference between fi u and fid is 
required. In the present case this is about 35 MeV. Pion 
condensation (i) and the 2SC phase (q) are in this dia- 
gram separated by two phase transitions in contrast to 
the estimated (fiB,fJ-i) phase diagram of Ref. |50| . 
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B. 



Hd = and T = 



In Fig. [2] we display the phase diagram for = 
and T = as a function of /x u and /x s . Since the up and 
down quark masses are much smaller than the strange 
quark mass, this diagram is very different from Fig. ^ 
Besides the possibility of pion condensation in (h) and 
(i), phases in which the charged kaon (k) and the neutral 
kaon condense (l)/(m) arise. The lines separating the 
charged kaon phase (k) from the chirally broken phase 
(a) are diagonal because at T — kaon condensation can 
occur if \fi s — fi u ,d\ > mK — 450 MeV [2a| (the chosen 
parameter set gives rise to a somewhat low kaon mass, 
but this is not relevant for the qualitative features of the 
phase diagram). In (r) we find the 2SCus phase. This 
phase is surrounded by phases in which the pions (h)/(i) 
and the neutral kaons (1) condense. We find that one 
passes a first-order transition when going from the pion 
and neutral kaon condensed to the 2SCus phase. 



H. (MeV) 




-600 -300 b 300 

\x u (MeV) 



600 



FIG. 2: Phase diagram for fid — and T — as a function of 
fi u and fi 3 . First and second-order transitions are indicated 
by solid and dotted lines, respectively. The letters denote the 
different phases, where a: uu + dd + ss, b: uu + dd, d: dd + 
ss, h: 7T+/7T-, i: tt+ /tt" + ss, k: K+/K~ + dd, 1: K°/K°, 
m: K°/K° + uu and r: 2SCus + dd. 



l^. fl u ~ fid 

In Fig. |3| we display the phase diagram for T = as 
a function of the up and down quark chemical potential 
and the strange quark chemical potential. We have cho- 
sen fi u = fid + e where e is a very small positive number. 
This e is necessary because when e = one is just at a 
first-order phase boundary between the phase in which 
the charged kaons condense (k) and the one in which the 
neutral kaons condense (m), as can be seen from Fig. EI 
This nonzero value of e gives rise to a small asymmetry in 



the phase diagram. If e is chosen negative, the phases in 
which the neutral (l)/(m) and the charged (j)/(k) kaon 
condenses are interchanged. 



600 



300 



H, (MeV) 



-300 




1 



-600 

-600 n -300 300 600 

Hu,Hd (MeV) 

FIG. 3: Phase diagram for T = as a function of fi u = fid + e 
and jj, s - First and second-order transitions are indicated by 
solid and dotted lines, respectively. The letters denote the 
different phases, where a: uu + dd + ss, b: uu + dd, j: 
K + /K~, k: K+/K~ + dd, 1: K°/K°, m: K°/K° + uu, n: 
2SC, q: 2SC + ss and t: CFL. 

Apart from the additional 2SC (n)/(q) and CFL (t) 
phases, our results agree qualitatively with the three- 
flavor calculations of Ref . [3^ (see their Fig. 7). 

The authors of Ref. y^] used different quark masses 
and a different coupling constant, and in addition em- 
ployed a form factor to mimic asymptotic freedom. 
Therefore, one may conclude that the use of such a 
form factor does not affect the phase diagram qualita- 
tively. We would also like to point out that the phase 
diagram Fig. [3] cannot simply be obtained by a super- 
position of phase diagrams obtained from a calculation 
with pseudoscalar condensation, but without supercon- 
ductivity (such as done in |32jp. and one with supercon- 
ductivity, but without pseudoscalar condensation (such 
as done in 51]). Despite the fact that the two types of 
phases do not coexist, there is nevertheless competition 
between them. Figure 03 shows that the K°/K° (l)/(m) 
and the K + /K~ (j)/(k) phases are separated from the 
2SC phase (q) by a first-order transition. This remains 
the case at finite temperature as is illustrated in Fig. 0] 
This figure displays the phase diagram as a function of 
fj, s and temperature, for fixed fi u = fid = 550 MeV. 

Returning to the discussion of Fig. [31 the line fi u = 
Md = Ms goes through the phase (a) in which chiral sym- 
metry is spontaneously broken. At some point it enters 
via a first-order transition the 2SC + ss phase (q), and 
finally goes into the CFL phase (t), again via a first-order 
transition. If there is a difference between fi u = /id and 
fi s , one can see from Fig.|3|that as the densities increase, 
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200 




-600 -300 300 600 

lis (MeV) 

FIG. 4: Phase diagram as a function of fi a and T, for fixed 
Hu = fid = 550 MeV. First and second-order transitions are 
indicated by solid and dotted lines respectively, while cross- 
overs are indicated by dashed-dotted lines. The uppermost 
phase without label is the restored phase. The letters denote 
the different phases, where g: ss, 1: K°/K°, n: 2SC, q: 2SC 
+ ss and t: CFL. 



quark matter can go directly from a phase of spontaneous 
chiral symmetry breaking (a) to a CFL phase (t) with- 
out passing the 2SC phase (q) first. This can also occur 
in compact stars |13|. One should keep in mind though 
that the relation between chemical potential and num- 
ber density is not linear. For example, at a first-order 
phase boundary, the number density increases discontin- 
uously. At these particular densities, quark matter can 
be in a mixed state of normal and superconducting mat- 
ter mm. 

It is also interesting to note that the phases (l)/(j) 
of kaon condensation can also occur outside the region 
of spontaneous chiral symmetry breaking. Assuming the 
phase transition towards chiral symmetry restoration co- 
incides with the deconfinement transition (as appears to 
be the case in lattice studies at small baryon chemical 
potential and in some models), this would imply that 
condensation of a state with quantum numbers of the 
kaon may persist in the deconfined phase. This was first 
observed in Ref. |27j , based on a perturbative calculation 
at high isospin chemical potential, that is expected to be 
applicable only in the deconfined region. 



D. fj, u = fi 3 

In Fig. \5\ we show the phase diagram at zero tem- 
perature as a function of fi u = fi s and fid- This dia- 
gram is similar to Fig. ^ for small strange quark chemi- 
cal potentials (below the kaon mass). At larger strange 
quark chemical potentials the diagrams differ, exhibiting 
kaon condensation (1) and diquark condensation involv- 



ing strange quarks (r)/(o)/(t). 




fj, u ,(i a (MeV) 

FIG. 5: Phase diagram for T — as a function of fi u = fi s 
and fid- First and second-order transitions are indicated by 
solid and dotted lines, respectively. The letters denote the 
different phases, where a: uu + dd + ss, c: uu + ss, d: dd + 
ss, g: ss, i: tt+ /tT + ss , 1: K°/K°, o: 2SCus, q: 2SC + ss, 
r: 2SCus + dd and t: CFL. 

In Fig. El we show the phase diagram as a function 
of fid and T, for fixed fi u = fi s = 550 MeV. In this 
figure one can find five critical points. Also, one can 
see in this figure that the K° /K° phase (1) is separated 
from the 2SCus phase (r) by a first-order transition for 
all temperatures. Furthermore, it is interesting that at 
finite temperature there is a first-order transition from 
the phase in which the neutral kaons condense (1) to the 
pion condensed phase (h). In Fig.0 we have enlarged the 
lower-right corner of Fig. [S] for clarity. In this figure one 
can find all possible superconducting phases, including 
the more exotic uSC (u), dSC (v) and sSC (w) phases. 
For fi u = fJ-d = fJ-s one goes from the CFL phase (t) via 
the 2SC phase (n) to the chirally restored phase when 
raising the temperature. However, small differences be- 
tween fi u — fi s and fid can cause one to go through com- 
pletely different phases. 



V. CONCLUSIONS 

In this paper, we studied the phase diagram of the 
three-flavor NJL model including pseudoscalar condensa- 
tion and color superconductivity as a function of the dif- 
ferent quark chemical potentials and temperature. The 
NJL model has a rich and interesting phase structure. 
The pseudoscalar condensed and color superconducting 
phases are competing and are separated by a first-order 
phase transition. As we have discussed, this need not 
be the case for other (less conventional) choices of the 
parameters of the model. 
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FIG. 6: Phase diagram as a function of /id and T, for fixed 
= Ms — 550 MeV. First and second-order transitions are 
indicated by solid and dotted lines, respectively, while cross- 
overs are denoted by dashed-dotted lines. The uppermost 
phase without label is the restored phase. The letters denote 
the different phases, where f: dd, h: tv + /tv~, 1: K°/K°, o: 
2SCus, r: 2SCus + dd, t: CFL and w: sSC. The lower right 
corner of this figure is enlarged in Fig. [7| 
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Furthermore, we concluded that at zero tempera- 
ture and zero strange quark chemical potential, there is 
a minimum asymmetry of about 35 MeV between the up 
and the down quark chemical potentials required in order 
to have two first-order transitions, when going from the 
phase with spontaneous chiral symmetry breaking to the 
2SC phase. 

Our results provide a qualitative check and extension 
of several earlier calculations that appeared in the liter- 
ature. The new aspects of the phase diagrams are often 
located in regions, where the quark chemical potentials 
are large and very different in magnitude for the differ- 
ent flavors. Although such situations are not necessarily 
realized in compact stars or can be realized in heavy-ion 
collisions, a comparison with future lattice data may nev- 
ertheless provide interesting information. This is espe- 
cially relevant for pseudoscalar condensation in the phase 
where chiral symmetry is restored and also for the com- 
plicated superconductivity phase structure close to the 
cutoff of the model. 

This work can be extended in several ways. For ex- 
ample, one can take into account 't Hooft's instanton- 
induced interaction |53[ . If one has pseudoscalar conden- 
sation, this is more difficult than in the normal case. An- 
other useful extension would be the inclusion of the neu- 
trality conditions 0, > m which case the phase struc- 
ture changes and for instance gapless phases will occur. 
It would also be interesting to see how the results depend 
on the strength of the diquark coupling and also on the 
quark masses. Furthermore, one could add the LOFF 
phase |45j. In this crystalline phase, quarks of different 
momenta can pair. One could also include vector inter- 
actions. In this case spin-1 diquark condensation (see for 
example Refs. 0>EU) an d an induced Lorentz-symmetry 
breaking (ISB) phase [5(| are among the possibilities. It 
would also be worthwhile to take pseudoscalar diquark 
condensation [34| into account. Finally, one could try to 
go beyond the mean-field approximation as was done in 
Ref. " 
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